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Abstract 

In this paper, we study the interplay between the epidemic spreading and the 
diffusion of awareness in multiplex networks. In the model, an infectious disease 
can spread in one network representing the paths of epidemic spreading (contact 
network), leading to the diffusion of awareness in the other network (information 
network), and then the diffusion of awareness will cause individuals to take social 
distances, which in turn affects the epidemic spreading. As for the diffusion of aware¬ 
ness, we assume that, on the one hand, individuals can be informed by other aware 
neighbors in information network, on the other hand, the susceptible individuals can 
be self-awareness induced by the infected neighbors in the contact networks (local 
information) or mass media (global information). Through Markov chain approach 
and numerical computations, we find that the density of infected individuals and 
the epidemic threshold can be affected by the structures of the two networks and 
the effective transmission rate of the awareness. However, we prove that though the 
introduction of the self-awareness can lower the density of infection, which cannot 
increase the epidemic threshold no matter of the local information or global informa¬ 
tion. Our finding is remarkably different to many previous results-local information 
based behavioral response can alter the epidemic threshold. 
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1 Introduction 


The outbreaks of diseases can involve the diffusion of information in regard 
to the diseases, including the risk of infection, rumors, fears and so on, which 
can stimulate individuals to take spontaneous behavioral responses to pro¬ 
tect themselves, thereby bring profound impacts on the spreading of dis¬ 
ease [T]|2]l3]ini5]l6] ■ For example, recent outbreaks of the HlNl flu, the bird 
flu, and the severe acute respiratory syndrome (SARS) have brought the 
reduction of going out and the plenty of people wearing face masks. For 
this reason, there has been an increasing focus on the development of for¬ 
mal models aimed at investigating the interplay of epidemic spreading and 
information-based behavioral responses [TIISlIMTO] . Such as, based on the as¬ 
sumption that the probability of susceptible individual going to the alter state 
is proportional to the number of infected neighbors, Sahnehef al. extended the 
SIS (Susceptible-Infected-Susceptible) model to a Susceptible-Alter-Infected- 
Susceptible (SAIS) model [TTIIT^ . and they found that the way of behavioral 
response can enhance the epidemic threshold; Melon! et al. studied a meta¬ 
population model that incorporates several scenarios of self-initiated behav¬ 
ioral changes into the mobility patterns of individuals, and they found that 
such behavioral changes do not alter the epidemic threshold, but may produce 
a negative impact on disease, i.e., the density of infection is increased [13]; In 
Refs. [milSIllG] . authors investigated the effects of the information-based be¬ 
havioral responses on the epidemic dynamics by designing the transmission 
rate as a function of the local infected density or the global infected density. 

Though the effects of information-based behavioral responses on the epidemic 
dynamics have been studied by many authors, most of works assumed the 
spreadings of information and epidemic are in the same network. As we know, 
with the development of technology, information can fast diffuse through many 
different channels, such as, the word of mouth, news media, online social net¬ 
works, and so on. In view of this, recent well-studied multiplex network theory 
has been used to mimic the interplay of information or related awareness and 
the epidemic dvnamics |17lll8lll9ll20] . For instance, Sahneh et al. have shown 
that the information dissemination spread in another network can help boost 
the resilience of the agents’ population against the spreading and found opti¬ 
mal information dissemination for different topologies [2l]; Wang et al. inves¬ 
tigated the interplay of the epidemic dynamics and the information dynamics 
in multiplex network based on the SIR (R-recovery) model, and focused on 
the two fundamental quantities underlying any spreading process: epidemic 
threshold and the hnal epidemic prevalence [ 22 ]; Granell et al. established 
an SIS-UAU model to investigate the competing effects of the spreading of 
awareness and the epidemic dynamics in multiplex with the transmission rate 
of awareness as well as the structure of information network [22]. More re¬ 
cently, they further generalized their model by reducing the probability of 
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infected individuals becoming awareness and including the effect of a mass 
broadcast of awareness (mass media) on the epidemic dynamics [23]. 

Existing works either assume that individuals are self-initiated uma, that 
is, they become aware because their neighbors are infected, or individuals can 
only be informed and become aware by other aware neighbors [221123] , the com¬ 
bine effects of the two factors have not been well studied. In reality, individuals 
can obtain the disease information through many ways. Inspired these factors, 
in the current work, we study the interplay between the diffusion of aware¬ 
ness by incorporating the self-awareness effects and the epidemic dynamics 
under the framework of multiplex network. In the model, an infectious disease 
hrst spreads among population represented by the contact network, and then 
the outbreak of the disease stimulates some people (infected or susceptible 
individuals) become aware of the risk of infection, and they take some protec¬ 
tions to reduce the probability of infection. Meanwhile, unaware individuals 
can be informed by other aware individuals through the information network 
or become self-awareness induced by the infected neighbors in contact net¬ 
work or mass media. The Ending indicates that the additional self-initiated 
awareness mechanism can reduce the density of infection, however, which can 
not alter the epidemic threshold. The results are verified by the Monte-Carlo 
simulations and the microscopic Markov chain approach (MMCA). 

The layout of the paper is as follows: we introduce the model in Sec. |2l The 
simulation results and theoretical analysis are presented in Sec. [3l Finally, 
Conclusions and discussions are presented in Sec. 01 The results for the global 
information-based self-initiated awareness are given in Appendix. 


2 Model 


In this work, we generalize the model of Ref. [23] . In that model, a multiplex 
includes two layers, one is physical layer representing the spreading of epidemic 
(contact network), and the other is information layer on where the diffusion 
of the awareness evolves (information network). All nodes represent the same 
individuals in both layers, but the connectivity is different in each of them. In 
the contact layer, a Susceptible - Infected - Susceptible (SIS) model is used 
to mimic the epidemic dynamics. That is to say, a susceptible node can be 
infected by one infected neighbor with certain probability, and the infected 
node can return to susceptible state with probability /x. On the information 
layer, the dynamical process of awareness is assumed to be similar to the SIS 
model, that is, an unaware node (U) can be informed by an aware neighbor (A) 
with probability A, and the aware node can loss awareness and back to unaware 
state with probability S. The interplay of the two processes is modelled as 
follows: once an individual is infected, s/he will certainly become aware, that 
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is, the probability is cr = 100%. In addition, to distinguish the protective 
behaviors between the aware individuals and unaware individuals, let {3 and 
{3^ = 7/5 (here 0 < 7 < 1 . If 7 = 0, the aware individuals are completely 
immune to the infection.) be the probabilities of unaware and aware susceptible 
nodes to get infected, respectively. 

From the description of the model, one can hnd that, on the one hand, the au¬ 
thors assumed that the infected individuals will automatically become aware 
and are willing to inform the disease information. As we know, in many cases, 
infected individuals are unwilling to tell others since they can be discrimi¬ 
nated or isolated by others once others know they are infected by one certain 
disease. So we make a progress and assume that infected individuals becom¬ 
ing aware with probability 0 < a < 1. On the other hand, in the model, 
individuals can only be informed by their neighbors through the information 
network, that is to say, one individual has no chance to become aware once 
their neighbors are unaware. However, as proposed in many previous works, 
individuals can become self-awareness once their friends are infected or they 
are informed by the mass media. Thus, we further generalize the model as: 
a susceptible individual can go to aware state by self-initiated response with 
probability k when contacting one infected friend. Therefore, the probability 
of the susceptible individual becoming awareness increases with the number 
of infected neighbors in the contact network [T2]. Note that, for the original 
model in Ref. [23], awareness cannot break out if the probability a = 0, so 
the roles of awareness cannot be played, however, in our model, the awareness 
can diffuse among population even though cr = 0 since susceptible individuals 
can become awareness by their self-initiated responses. When k = 0.0 and 
cr = 1.0, our model returns to the original model in Ref. [23] . 

According to this scheme, an individual can be in four different states: suscep¬ 
tible and unaware(SU), susceptible and aware(SA),infected and unaware(IU), 
infected and aware(IA). The flow diagram of the model is given in Fig. [H 


3 Main results 


3.1 Theoretical analysis 


Denoting Ojj and bij be the adjacency matrices that support the SIS and 
UAU processes, respectively. The probability of i in one of four states at time 
t is denoted by pf^(t), pf^(t), pl^ {t) and pl^{t) respectively. Assuming the 
probability of susceptible (infect) node i not being informed by any neighbors 
is 9i(t) (rj(f)), and the probability of unaware (aware) susceptible node i not 
being infected is qY(t) {qf(t)). They are described as [23] : 
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Fig. 1. Model description for the UAU-SIS dynamic. An individual can be in four 
different states; SU, SA, lU, and lA. The top (bottom) layer is the epidemic process 
for the aware (unaware) individuals, respectively. SU (SA) can be infected by an 
infectious neighbor in contact layer with a probability /3 (/3"^ = 7/3). lU and lA re¬ 
covers to SU and SA, respectively, with the same probability /r; The left (right) flow 
is the awareness process for the susceptible (infected) individuals, respectively. SU 
can go to SA with a probability A of being informed by an aware neighbor through 
information network, or induced by the infected neighbors in contact network with 
a probability k. SA recovers to SU with a probability 6. lU can go to lA by informed 
the aware neighbors in information layer with a probability A, or self-awareness with 
a probability a. SA and IA can become unaware and return to SU and lU with the 


! probability 5. 


=11(1 - bjipf{t)X){l - ajip]{t)K), 

j 

(1) 

ri{-t) =Y[i^ - bjipf{t)\), 

j 

(2) 

j 

( 3 ) 

j 

( 4 ) 


where pf(t) = Pj^(t) + Pj^(t) and Pj(t) = Pj^(t) + Pj^it) (Note: to simplify 
the model, we do not distinguish the infectivity of lA and lU, meanwhile, the 
diffusion capabilities of awareness for SA and lA are also the same.). 
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Fig. 2. Transition probability trees for the states of UAU-SIS dynamics in the multi¬ 
plex for per time step. The denotations of 9i, ri, and qf are given in Eqs. m-m- 
Other parameters have the same denotations as in Fig. [1] 

For each possible state at time t, which may give rise to four possible states 
at time t + 1 with certain probability, the transition probability trees for node 
i are illustrated in Fig. |2l According to the scheme in Fig. [2l one can easily 
write the Markov Chain Approach (MMCA) equations [2^26] for each state, 
for example, along the top branch of the four trees, one can read the probability 
pf^{t + 1) of SU at f + 1 as: 


pf^{t + 1) = pl^{t)6p + pI'^ + pf^{t)5q'^ (t) + pf'^ {t)9i{t)q'/ (t)(5) 


JA/ 


JU 


SA/ 


JU/ 


M/ 


Similarly, the other three MMCA equations can be written as: 


+ 1 ) - ri{t))p ( 6 ) 

+ pf 

+ 1) = (1 -/i) -/i) (7) 
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+pf^{t)6{l - - gf (t))|, 

+ 1) = - p)+pl^it)ri{t){l - p)+pf^{t)6{l - gf (t)) 

+pf^{tWi{t){i - q^{t))^ +pl^{t){l - 5)(1 - p) (8) 

+pi^{t){l - ri{t)){l -p)+ pf^{t){l - 5)(1 - qf{t)) 

+Pf^it)i^ - Oiit))il - qf{t)), 

where pf^{t) +Pi^{t) +Pi'^{t) +Pi^{i) = 1- When the system goes to the sta¬ 
tionary state, we have pf^{t + l) = pf^it) = pf^ for SU state and equivalently 
for the SU, SA and lU states. 

Since the epidemic threshold determines whether the epidemic can outbreak 
or die out, it is vital to analyze the effects of the different parameters on 
the epidemic threshold /5c- As stated in Ref. [23], near the threshold, the 
probability of nodes being infected is very low, i.e., pi = £i ^ 1. Consequently, 
~ 1 — (3^ Y^jicLjiSj) and qY ~ 1 — /5 Y^j^ajiSj). Further approximate qf ~ 1 
and gf 1 by assuming pi = Si ^ 0, and then substitute both of them into 
Eq. (|5|) and Eq. (|6|), we obtain 

pf = pf^s + pfe, (9) 

py=?“(!-^)+pf(i-Si) (10) 


Combing Eqs. (7-10), then a simple formula is obtained: 

= {pf^P + Il(aiUi)- (11) 

With (3^ = 7/5, pY = + pI^ ~ pf^, pf = pf^ + pl^ ~ pf^ and Ei = 

Pj’^it) +Pj'^it) <C 1, then Eq. flTTl) can be rewritten as: 

E[(l - (1 - 7)pf )ap - = 0, (12) 

here 6ji = 1 if i = j; otherwise, dij = 0. 

Dehning matrix H with elements: 

hji = (!-(!- j)pf)aji, (13) 

Eq. (fT2]) can be read as 

He = (14) 
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where £ = (ei, £ 2 , • ■ • ,£n)'^ with T be the vector transportation. 


The non-trivial solntions of Eq. flTT|) are eigenvectors of H, whose eigenvalnes 
are eqnal to /r//9. Therefore, the onset of the epidemics is given by the largest 
eigenvalne of A^axiH), 


Ac 


h 

Amax ('^) 


(15) 


3.2 Numerical simulations 


To verify onr theoretical resnlts, as in Ref. [23], we bnild a confignration 
network with degree distribntion P{k) ~ and network size N = 2000 as 
the contact network, and for the information network, which is generated by 
adding 800 extra random links in the contact network, and Aj represent 
the adjacency matrices of the contact network and the information network, 
respectively. All simnlation resnlts are obtained by averaging 20 realizations. 

We first compare the resnlts from MMCA with Monte-Carlo simnlation in 
Fig. [3] to check the effectiveness of onr analysis based on MMCA, from Fig. [3l 
one can observe that the resnlts based on the two approaches are in good 
agreement. So in the next hgnres, onr main results are obtained from MMCA. 

Then we investigate the effects of the two main parameters of the model— 
K and a on the epidemic threshold and the density of infected individuals. 
Here, we will present the results for 7 = 0, meaning that = 0 and q'^ = 0. 
Obviously, once the value of 7 is increased, the epidemic threshold is decreased 
and the density of infected individuals in enhanced correspondingly. 

Fig. m plots the density of infection [pA see Fig. 111(a)] and aware individuals 
[p^, see Fig. 111(b)] as a function of A for different values of k, respectively. 
Observing the figure, one can see that though the larger value of k, can cause 
more individuals become aware and reduce the density of infection. However, 
one can find that the increasing of A has no influence on the epidemic thresh¬ 
old. The result is remarkably different from many previous results which claim 
that the local information-based behavioral response in the single-layer net¬ 
work can alter the epidemic threshold. How to understand this nontrivial 
result? Since the UAU awareness dynamic is the same to the SIS epidemic 
process. When the epidemic has not broken out, the density of awareness in 
information network (he., p"^) is only determined by the transmission rate of 
awareness. A, recovery rate 6 and the structure of information network, but is 
not related to the value of k. Namely, near or below the epidemic threshold 
point, increasing the value of k, only means that the initial number of the aware 
individuals is increased, which cannot affect the density of aware individuals 





Fig. 3. Comparison of MMCA(solid line) with Monte Carlo simulation (dotted line) 
for A = 0.15, 7 = 0, 5 = 0.6 and [i = 0.4. The fraction of infected (a) and aware 
(b) nodes as a function of the infectivity parameter (3 for two different conditions 
of the parameters k and a. 

at stationary state. In this case, the valne of pf is independent of the value of 
K or a, which gives rise to the same value of [see Fig. 111(b)]. Thus, accord¬ 
ing to Eqs. flT^ and flT^ . the epidemic threshold f3c is invariable owing to the 
same value of p^. One should note that, for the case of single-layer network, 
the local information-based behavioral response can directly reduce the trans¬ 
mission rate of epidemic, leading to the change of the epidemic threshold. For 
our model, the self-awareness behavior Erst diffuses through the information 
network, and then the epidemic process happens in contact network. Thus, 
the effect of the awareness behavior on the transmission rate is indirect. 

The density of infected individuals and aware nodes as functions of (3 for 
different values of parameter a are also shown in Fig. |5](a) and Fig. l5](b), 
respectively. One can see that, similar to Fig. HI varying the value of a has 
no effect on the epidemic threshold. However, differ to the above case, from 
Fig. [5](a) we find that the value of a also has negligible effect on the density of 
infection, even in the extreme cases where infected unaware individuals remain 
unaware of its sickness {a = 0) or certainly become aware of it {c = 1). The 
result in Fig. [5] is in accord with the Fig. 3 in Ref. [23]. 

In order to systematically study the effects of k and a on the we further 
explore the full phase diagram (A — (3) in Fig. [6l Overall, we can see that p^ 
is not influenced by A when f3 is smaller than the epidemic threshold, since 
epidemic will die out by itself. Once (3 overpasses the epidemic threshold, p^ 
decreases with A for different values of k or a. More specifically, by comparing 
Fig. |6](a) with Fig.[6](b) (or comparing Fig. |6](c) with Fig. |6](d)), we can see that 
p^ is not remarkably influenced by the value of a. Likewise, by comparing the 
Fig. |6](a) with Fig. |6](d) (or comparing Fig. |6](b) with Fig. [6](c)), p^ decreases 
with K, especially for the large value of (3. 
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Fig. 4. Fraction of infected (a) and aware (b) nodes as a function of the infectivity 
parameter /3 for different values of the parameter k. Here A = 0.15, 7 = 0, <5 = 0.6, 
/i = 0.4 and a = 0.4. 




Fig. 5. Fraction of infected (a) and aware (b) nodes as a function of the infectivity 
parameter /3 for different values of the parameter a. Here A = 0.15, 7 = 0, 5 = 0.6, 
// = 0.4 and k = 0.05. 

4 Conclusions and discussions 


Recognizing that, on the one hand, ontbreak of an epidemic throngh a physical- 
contact network can trigger the spreading of information awareness through 
other different channels, such as on-line social networks, mass media; on the 
other hand, an individual can not only be informed by other aware neighbors 
but also can become self-awareness once some friends in contact network are 
infected. By introducing the self-awareness mechanism for susceptible individ¬ 
uals, we have investigated interplay between the spreading of epidemic and the 
diffusion of awareness based on the framework of the multiplex networks. We 
mainly studied the two parameters k and a characterizing the self-awareness 
probability of susceptible individuals and infected individuals, respectively. 
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Fig. 6 . The fraction of infected individuals in the stationary state. Full phase 
diagram A — /3 for the same multiplex described before Where 7 = 0, A = 0.15, 
5 = 0.6 and p = 0.4. (a): k = 0, a = 1; (b): k = 0, a = 0.2; (c): k = 0.2, a = 0.2; 
(d): K = 0 . 2 , a = 1. 

Analysis based on the Markov chain approach as well as the extensive com¬ 
putations reveal that the density of infection can be reduced once the two 
parameters are increased, however, we found that the impact of self-awareness 
behavior for susceptible individuals on inhibiting the spreading of epidemic 
is much better than the self-awareness of the infected individuals, since self- 
awareness from susceptible individuals can directly reduce their probabilities 
of being infected. Meanwhile, we found that the self-awareness behavior can¬ 
not alter the epidemic threshold no matter of the local or global information, 
which are in stark contrast with the results obtained from the single-layer 
networks. 

The challenges of studying the intricate interplay between social and biological 
contagions in human populations are generating interesting science. In this 
work, we consider the effects of the self-awareness behavior based on the 
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multiplex networks on the density of infection and the epidemic threshold, 
our result implies that the conclusions obtained from single-layer networks 
may need to be re-examined when they are extended to multiplex networks. 
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6 Apendix: Global information-based awareness 


In recent work, Granell et al. considered the effect of the mass media on the 
epidemic process and awareness diffusion |23]. In the model, each individual 
becomes aware with probability m by assuming that they are informed by a 
broadcast or mass media. Thus, it can be regarded as a global information- 
based awareness. One questionable point is that the probability of being aware¬ 
ness m is irrelevant to the density of infection. As a result, even the epidemic 
is almost eliminated, individuals also have the fix probability of being aware. 
In reality, becoming awareness often means that individuals need to take some 
protective measures, such as, washing hands, wearing masks or reducing out¬ 
goings. These measures indicate certain inconveniences or some cost [2Z12B]. 
Thus, a more realistic situation is that the probability of being awareness 
should be related to the density of infection. To mimic this case, here we as¬ 
sume the probability of aware from global information is given as: ml[t) with 
I{t) is the density of infection at time t, which indicates that the probability 
of being awareness adaptively varies with the density of infection. 

For this case, we only need to slightly change the local model described in 
subsection 13.11 We only need to change 'jiit) and 9i(t) as follows: 

= n(l “ bjipf{t)X){l - m/(f)), (16) 

j 

= n(l “ - mp\t)) (17) 

j 


Similarly to the above analysis, we can get that the epidemic threshold is still 
determined by Eq. (fT5|) . i.e., the epidemic threshold is also independent of the 
value of m, which is different from the result in Ref. [23] . The result is verihed 
by Fig. [7] for different values of 7 . 
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Fig. 7. The density of infected individuals as a function of f3 for different values of 
m. Here A = 0.3, a = 0.5 6 = 0.6, ^ = 0.4 and k = 0.0. (a): 7 = 1; (b): 7 = 0.5; (c): 
7 = 0.25; (d): 7 = 0 . 
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